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ABSTRACT 

Aims. We study the close circumstellar environment of the nearby S-type star -k 1 Gruis using high spatial-resolution, 
mid-infrared observations from the ESO/VLTI. 

Methods. Spectra and visibilities were obtained with the MIDI interferometer on the VLT Auxiliary Telescopes. The 
cool M5III giant j3 Gruis was used as bright primary calibrator, and a dedicated spectro-interferometric study was 
undertaken to determine its angular diameter accurately. The MIDI measurements were fitted with the ID numerical 
radiative transfer code DUSTY to determine the dust shell parameters of -k 1 Gruis. Taking into account the low spatial 
extension of the model in the 8-9 (im spectral band for the smallest projected baselines, we consider the possibility of 
a supplementary molecular shell. 

Results. The MIDI visibility and phase data are mostly dominated by the spherical 21 mas (694 Rq) central star, while 
the extended dusty environment is over-resolved even with the shortest baselines. No obvious departure from spherical 
symmetry is found on the milliarcsecond scale. The spectro-interferometric observations are well-fitted by an optically 
thin (rdu S t<0.01 in the N band) dust shell that is located at about 14 stellar radii with a typical temperature of 700 K 
and composed of 70% silicate and 30% of amorphous alumina grains. An optically thin (r mo i<0.1 in the N band) 
H20+SiO molecular shell extending from the photosphere of the star up to 4.4 stellar radii with a typical temperature 
of 1000 K is added to the model to improve the fit in the 8-9 fim spectral band. We discuss the probable binary origin 
of asymmetries as revealed by millimetric observations. 

Key words. Techniques: interferometric; Techniques: high angular resolution; Stars: AGB and post-AGB; Stars: circum- 
stellar matter; Stars: mass-loss 



\ 1. Introduction et al., 1998), tt 1 Gruis (Hen 4-202, HD212087, HIP 110478) 

■ . . r n . . . . _ _ , . is one of the brightest intrinsic S stars. It is a semi-regular 

"£j . The class c,f S stars is characterized by ZrO (and some- variable Qf gRb . from y=5 4 tQ g ? ^ 15Q 

; times LaO) bands, m addition to the TiO bands present d It wag chosen ag an Qri ^ tot of the s star 

m cool oxygen-rich giants of spectral type M. The presence ^ ^ R And and R c (Men . m m2) and ig knQwn 

ot lines from the unstable clement 1c (technetium m the , , u ■ , , , r<r\\r ■ i 

. . . . _ ... v ' to be a binary star with a solar-type GOV companion, lo- 

spectra of intrinsic S stars indicates that they are currently cated ftt &n a kf digtance Qf 2 n „ (unch a d for more 

evolving on the thermally-pulsmg asymptotic giant granch . r> ■ \ <a.t. -j.- 

A ^„ x , . J , , ,i . . , than a century; Bonneau, priv. comm.), with a position 

TP-AGB), where internal nucleosynthesis is responsible k of 2Q2 r and an rent yisual magnitude of 10 . 9 

for forming carbon and heavy e ements (s- elements) sub- (Feag 1953) _ The k ig obabl h ical with their 

sequently dredged up to the stellar surlacc. Being on the , '., , • . . , . t 

. !1_ . J . . °„ * . , ,, . , parent magnitude agreeing with their spectral type tor the 

AGB, intrinsic S stars experience a mass- loss allowing the tt- a- ± 

r ' r i ■ i , • ii • , /t ■ Hipparcos distance, 

formation ot tlieir dust circumstellar environment (Jorissen 

& Knapp, 1998). With c / ratios generally around 0.8 (Smith & 

At a distance of 153+ 20 pc (Perryman et al., 1997; Van Eck Lambert, 1986), photospheres of S-type stars are oxygen- 

~ ! ~ ~~ \ \ „ „ rich. The situation regarding tt 1 Gruis is somewhat 

bend offprint requests to: S. Sacuto r . , A i- , • • TT ^ i 

* „ , , ,. , ... ,, ,. . _ . contusing, though. According to its position m IK color- 

Based on observations made with the Very Large telescope . ,. i „ , . , .„,. , „ , n 

Interferometer at Paranal Observatory under programs 077.D- color diagrams, tt Gruis is classified m the C class of 

0294(D/E/F) intrinsic S stars (Jorissen & Knapp, 1998). This class 

Reduced visibilities and differential phases are avail- includes stars with a dense oxygen-rich circumstellar enve- 

able in electronic form at the CDS via anonymous !ope and with strong chemical peculiarities and mass-loss 

ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via rates ranging from a few 10~ 7 to a few 10~ 5 Moyr^ 1 . 

http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/?/? Furthermore, this star is classified as E (emission silicate) 
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by Volk & Cohen (1989) and SE (silicate and oxygen-rich 
dust emission) by Sloan & Price (1998). The circumstellar 
dusty shell of n 1 Gruis would thus appear to be mainly 
composed of oxygen-rich elements; however, the track 
in the (K- [12], [25]- [60]) color-color diagram during the 
variability cycle (Jorissen & Knapp, 1998) suggests that 
the shell emission could as well be associated with carbon- 
composite material (Thomas et al., 1976; Knapp et al., 
1999a). 

The mass-loss history of the star can be traced back 
21,000 years, with the detection of a 0.28 pc extended dust 
shell by IRAS, assuming a 15 kms" 1 velocity of expansion 
(Young et al., 1993). The molecular environment was re- 
solved in the east-west direction with a 24" (Sahai, 1992) 
and a 30" beam (Knapp et al., 1999a). Moreover, the CO 
lines exhibit an unusual asymmetric double-horned profile 
and a high-velocity component interpreted in the frame of 
a slowly expanding disk and a fast (v > 90 kms -1 ) north- 
south bipolar outflow. Chiu et al. (2006) report millimeter 
observations at higher spatial resolution (2.1" x4. 2" beam) 
that confirm this interpretation. In addition, these observa- 
tions reveal a 200 AU cavity suggesting that the mass-loss 
rate of ir 1 Gruis has dramatically decreased over the past 
90 years. Whether the highly structured far environment 
of 7T 1 Gruis (presence of disk + jet) is a consequence of 
the asymmetry of the dusty wind from the primary star 
or of the presence of the far companion, or even a closer 
undetected one, cannot be addressed by the millimetric ob- 
servations of Chiu et al. (2006) due to insufficient spatial 
resolution. 

Based on the comparative analysis of the jets and disks 
present in a sample of late AGB stars (among which V Hya 
and 7T 1 Gruis) and of proto-planetary nebulae (including 
symbiotic stars), Huggins (2007) strongly favors the 
second possibility. The properties of the disk and bipolar 
outflow of 7T 1 Gruis appear to fit nicely in a sequence 
going from late (binary) AGB stars to proto-planetary 
nebulae. According to the binary scenarios discussed by 
Huggins (2007), the GOV companion could be at the 
origin of the launching of the bipolar outflow, despite 
the large separation. However, Makarov & Kaplan (2005) 
and Frankowski et al. (2007) found possible evidence for 
yet another companion from the discrepancy between the 
Hipparcos and Tycho-2 proper motion (the Hipparcos 
proper motion relics on measurements spanning only 
3 years, whereas the Tycho-2 one uses positional data 
spanning more than a century). Any discrepancy between 
these short-term and long-term proper motions is thus in- 
dicative of the presence of an unrecognized orbital motion 
in addition to the true proper motion. The discrepancy 
is very significant in the case of tt 1 Gruis. According to 
the analysis of Frankowski et al. (2007), these so-called 
'A/z binaries' correspond to systems with orbital periods 
in the range 1500-10 days, much shorter thus than the 
6000-year period inferred by Knapp et al. (1999a) for the 
GOV companion 2.7" away. It is therefore an open question 
whether a third component might be present in the system. 

In this paper, we present the first mid-infrared interfero- 
metric observations of this star taken with the MIDI/VLTI 
instrument, secured with auxiliary telescopes (ATs) of 
1.8 m that provide a beam with an FWHM of 1.1" at 
10 fim. In interferometric mode, measurements of spatially 



correlated flux at a high spatial resolution (^20 mas) and 
multi-spectral (A/AA~30) information on the mid- infrared 
source was obtained within the interferometric field-of-view 
(FOV) of 1.1". The inner diameter of the millimetric disk 
(^2.6") found by Chiu et al. (2006) is then largely beyond 
the present MIDI FOV, as is the distant companion (2.71"). 
The outline of this paper is as follows. In Sect. 2, we 
present the MIDI observations and describe the calibration 
performed with the two data-reduction software packages 
(MIA and EWS) developed for the MIDI focal instrument. 
This calibration required considerable effort to determine 
the angular diameter of the spatially resolved calibrator 
M5III star (3 Gruis (described in Appendix A). Section 3 
presents our spectro-photomctric and interferometric inter- 
pretations of the star. This was done first with an analytical 
radiative transfer model which gives some insight into the 
expected behavior of the mid-infrared emission from the 
circumstellar envelope. Then, we used the numerical radia- 
tive transfer code DUSTY to confirm the set of parameters 
found with the analytical thin-dusty-shell model and define 
a chemical composition of the dust shell, taking the posi- 
tion of the star on the silicate dust sequence into account. 
Moreover, to account for the larger than expected spatial 
extension of the object in the 8-9 /im spectral band, we 
add a molecular shell to the model. In Sect. 4, we discuss 
the probable binary origin of asymmetries revealed by mil- 
limetric observations. Finally, we summarize and conclude 
in Sect. 5. 

2. Observations and calibration process 

The Very Large Telescope Interferometer (VLTI) of ESO's 
Paranal Observatory was used with MIDI, the MID- 
infrared Interferometric recombiner (Leinert et al., 2003). 
MIDI combines the light of two telescopes and provides 
spectrally resolved visibilities in the N band atmospheric 
window. Guaranteed time observations (GTO) of n 1 Gruis 
were conducted with the VLT ATs positioned at stations 
AO, E0, DO, and GO, providing baselines in the range of 
15-64 meters. Observations were made during the nights 
of May 21-23-24-25-27 2006, June 19, 2006, and August 
8, 2006, under good atmospheric conditions (average 
seeing=0.94"), with the worst seeing (=1.51") during the 
night of May 27, 2006. Table 1 presents the journal of the 
interferometric observations. The calibrator (3 Gruis was 
either observed right before or after each science target 
observation. 

Chopped acquisition images were recorded (f=2Hz, 2000 
frames, 4 ms per frame, 98mas per pixel) for the fine acqui- 
sition of the target. The acquisition filter was an N-band 
filter. Simultaneous interferometric and photometric mea- 
surements were performed with the SCI-PHOT mode of 
MIDI, using the prism that provides a spectral dispersion 
of about 30. The data-reduction software packages 1 MIA 
and EWS (Jaffe, 2004) were used to reduce the spectra 
and visibilities (Chesncau et al., 2005). MIA is based on 
the power spectrum analysis and uses a fast Fourier trans- 
formation (FFT) to calculate the Fourier amplitude of the 
fringe packets, while EWS uses a shift-and-add algorithm 
in the complex plane, averaging appropriately modified in- 

1 http : //www . mpia-hd . mpg . de/MIDISOFT/ , 
http : / /www . strw . leidemmiv . nl/~nevec/MIDI/ 
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Table 1. Journal of observations. 



Star 


UT date &z Time 


4>* 


Config. 


B [m] 


PA [dc 


7T 1 Gruis 


2006-05-21 


09:09:26 




E0-G0 


15.7 


50 


(3 Gruis 


2006-05-21 


09:28:33 


0.23 


- 


- 


- 


7T 1 Gruis 


2006-05-23 


09:32:27 




A0-G0 


63.5 


56 


(3 Gruis 


2006-05-23 


10:06:51 


0.28 


- 


- 


- 


7r Gruis 


2006-05-24 


09:27:34 




D0-G0 


31.8 


56 


(3 Gruis 


2006-05-24 


09:44:39 


0.31 


- 


- 


- 


n 1 Gruis 


2006-05-25 


08:09:13 




A0-G0 


61.2 


41 


(3 Gruis 


2006-05-25 


07:51:28 


0.34 








7T 1 Gruis 


2006-05-27 


09:20:41 




E0-G0 


15.9 


57 


(3 Gruis 


2006-05-27 


09:37:11 


0.39 








7T 1 Gruis 


2006-06-19 


10:16:01 




D0-G0 


30.6 


83 


/3 Gruis 


2006-06-19 


09:50:14 


0.01 








7T 1 Gruis 


2006-08-08 


08:40:43 




A0-G0 


52.6 


102 


/3 Gruis 


2006-08-08 


08:13:56 


0.36 









* : the pulsation 
was not possible 
of its lightcurve. 



phase for j3 Gruis 
for -k 1 Gruis due 



(Otero & Moon, 2006). This 
to the non-regular character 



dividual exposures (dispersed channeled spectra) to obtain 
the complex visibility. 

An extensive study of the fluctuations of the instrumental 
transfer function was carried out for the baseline AO-GO 
during the night of 2006 May 25 (see Table A.l). The rms 
scatter of these calibration measurements is ^0.05 between 
8 nm and 13 /im, corresponding to a relative error of 10%. 
This error in SCI-PHOT mode with ATs is greater than the 
typical errors encountered with UTs with the same observ- 
ing mode (2-5%). This discrepancy is mostly due to the 
limited accuracy of the photometric extraction with ATs 
versus UTs. 

All visibility and differential phase data as well as all the 
characteristics of the observations are available from the 
CDS (Centre de Donnees astronomiques de Strasbourg); 
all data products are stored in the FITS-based, optical in- 
terferometry data exchange format (OI-FITS) described in 
Pauls et al. (2005). 

2.1. General description 

Due to its proximity (153 pc) and its brightness (about 900 
Jy at 10 /im), n 1 Gruis is an ideal mid-infrared target for in- 
terferometric measurements with the ATs. Interferometric 
observations provide the visibility measurements of the 
source V sou , which is the product of the source calibrated 
visibilities V sou with the system response V sys , which in turn 
is a combination of the atmospheric and instrumental re- 
sponses. The system response is estimated using the ratio 
between visibility measurements of the calibrator V ca i and 
its theoretical visibilities V ca \ (Van Belle & Van Belle, 2005). 
The expression of the source calibrated visibility is 



-v cal . 



(1) 



cal 



If the calibrator is unresolved, its theoretical visibility is 
close to 1 and the measured contrast of the fringes is a di- 
rect measure of the instrumental efficiency. If the calibrator 
is resolved, then the values of its theoretical spectrally dis- 
persed visibilities need to be very accurately determined. 
For our observations of ir 1 Gruis, which has no nearby un- 
resolved calibrator with enough flux (>100 Jy at 12 /mi, 



corresponding to the sensitivity of MIDI with ATs in low 
.spectral resolution and SCI-PHOT mode), two choices were 
_possible: either move the telescope far from tt 1 Gruis to 
acquire a bright enough unresolved calibration star or to 
observe a bright but resolved calibrator in the vicinity of 
7T 1 Gruis. The first choice introduces a bias due to the air 
mass difference between the source observation at a given 
sky location and the calibrator observation at another lo- 
cation. The second choice requires an accurate knowledge 
of the angular diameter of the resolved target. The second 
option was chosen for optimizing the flux at the price of 
a potentially large uncertainty on the angular diameter of 
the calibration star (see Appendix A). 

2.2. Calibrated flux of the source 



The MIDI flux of tt 1 Gruis was scaled to the MARCS model 
flux of (3 Gruis, for which we take into account the value 
of the diameter from the spectrometric method (26.8±1.3 
mas; see Appendix A.l) 



1 TTimarcs 
"^cal •> 



(2) 



cal 



where F sou is the calibrated flux of n 1 Gruis, F sou the raw 
flux measurements of tt 1 Gruis, F ca \ the raw flux measure- 
ments of (3 Gruis, and F™f rcs the MARCS synthetic flux of 
(3 Gruis. 

Figure 1 compares the single-dish MIDI spectrum with 
the ISO spectrum (Sloan et al., 2003) taken on 1996 Oct 
26. Error bars are related to the uncertainty on the value of 
the calibrator diameter. The MIDI flux level is lower than 
the ISO/SWS flux level. Two possibilities could account for 
this difference. The most probable one is that the beam of 
an AT (~1.1" at 10 /im) is small compared to that of the 
ISO/SWS satellite (14"x22"). Therefore, because of the 
proximity of this star (153 pc) and the large extension of 
its dust envelope beyond the AT beam size, some flux might 
be missed. A second possibility that cannot be excluded is 
that the level of emission has changed since the ISO ob- 
servation (1996) due to the pulsation of the star. However, 
this star is a semi-regular SRb pulsator with a low vari- 
ability, for which it is difficult to assign a phase to a given 
observation. Nevertheless, we can see that the shape of the 
ISO/SWS and MIDI spectra is about the same, revealing a 
dust feature that will be discussed in Sect. 3.4.2. As a con- 
sequence, we assume in the following study that the MIDI 
spectrum can be complemented with the ISO spectrum to 
construct the current SED of the object at longer infrared 
wavelengths. 

2.3. Calibrated visibilities of the source 

In Fig. 2 are shown simultaneously the calibrated visibility 
with the MIA and EWS reduction softwares. The curves 
from the two softwares are almost indistinguishable. For 
clarity, we only show the reduction taking the spectromet- 
ric diameter of (3 Gruis into account (26.8±1.3 mas; see 
Appendix A.l). The error bars are the sum of the trans- 
fer function error bar of 10% plus the error bar related 
to the uncertainty on the value of the calibrator diameter 
(from 0.6% at low spatial frequencies to 10% at high spatial 
frequencies). 
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Fig. 1. Calibrated MIDI fluxes of -k 1 Gruis (vertical error bars) 
as compared to the ISO/SWS spectrum (solid line) taken on 
1996 Oct 26. 
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indicating that an extended mid-infrared structure sur- 
rounds the central star (Ohnaka et al., 2005). The di- 
mension of this structure depends on the observed wave- 
length, so neither of the two simple models (monochro- 
matic uniform-disk or Gaussian distribution) can describe 
the source. In the following, we model the visibilities with 
radiative transfer models taking into account the spectral 
dependency of the circumstellar environment of tt 1 Gruis. 

To check whether the object shows some asymmetries, 
one can use the differential phase of tt 1 Gruis obtained from 
the EWS software. Figure 3 shows the calibrated differen- 
tial phase of the object evaluated from all the determined 
projected baselines. Error bars correspond to the rms of the 
differential phase over all the projected baselines. The very 
low rms values of the differential phase (±1.93° at 8 /zm, 
±1.10° at 9 fim, ±1.65° at 10 fim, ±1.82° at 12.5 /xm) in- 
dicates that there is no obvious signature of an asymmetric 
source, such as binarity or a clumpy environment, for the 
given projected baselines. 
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Fig. 3. Calibrated differential phase of n 1 Gruis evaluated from 
all the determined projected baselines (dotted lines). Error bars 
correspond to the rms of the differential phase averaged over all 
the projected baselines. 
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Fig. 2. MIDI-calibrated visibilities (error bars) with the MIA 
(squares) and EWS (triangles) reduction softwares. The solid 
lines correspond to the uniform disk angular radius (in mas) 
and the dashed lines to the FWHM (in mas) of the Gaussian 
distribution (both to be read from the scale on the right axis), 
computed from the visibility at each spectral channel. 



Furthermore, the MIDI visibility amplitudes also show 
no evident departure of the object from spherical symme- 
try. Figure 4 shows the dependence of the dimensions of 
the equivalent uniform-disk on the projected baseline po- 
sition angle on the sky at different wavelengths. The three 
circles fit the points calculated for a given wavelength at 
close projected baselines for different projected baseline po- 
sition angles. In addition to zero differential phase, these 
fits show that the object is nearly spherical. It must be 
pointed out that for such a close, extended giant, the in- 
terferometric fringes essentially probe the star photosphere 
and its close vicinity, while the extended dusty environment 
is over-resolved even for the shortest baselines. 



The size derived from the uniform disk and the Gaussian 
distribution regularly increases from 8 to 13 /*m (see Fig. 2), 



From these considerations, spherical radiative transfer 
models will be used to interpret the spcctro-intcrferometric 
data of tt 1 Gruis. 
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Table 2. Broadband photometry of it 1 Gruis 



-50 -40 -30 -20 



10 10 20 30 40 
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Fig. 4. This figure presents the dependence of the dimensions of 
the equivalent uniform-disk (in mas) on the projected baseline 
position angle on the sky taken at 8 (plus), 9 (cross), 10 (tri- 
angle), and 12.5 (im (square). The three dotted circles fit the 
points calculated for a given wavelength at close projected base- 
lines (B) for different projected baseline position angles (PA). 



Band 


Wavelength 


Extinction 


Dereddened 


Reference 




H 


coefficient 


magnitude 




R 


0.67 


0.13 


3.05 


(1) 


I 


0.87 


0.09 


0.45 


- 


H 


1.65 


0.03 


-1.82 


(2) 


K 


2.2 


0.02 


-2.07 


- 


L 


3.6 


0.01 


-2.50 


- 


M 


4.8 


0.01 


-2.26 


- 


N 


8.4 


0.01 


-2.89 


- 


- 


9.7 


0.01 


-3.25 


- 


- 


10.5 


0.01 


-3.51 


- 


- 


11.2 


0.01 


-3.58 


- 


- 


12.5 


0.0 


-3.47 


- 


Q 


20.0 


0.0 


-4.35 


- 


u 


0.35 


0.28 


11.13 


(3) 


B 


0.42 


0.24 


8.32 


- 


V 


0.55 


0.17 


6.43 


- 


J 


1.21 


0.05 


-0.60 


- 


H 


1.66 


0.03 


-1.72 




K 


2.22 


0.02 


-2.14 




L 


3.45 


0.01 


-2.57 




J 


1.25 


0.05 


-0.71 


(4) 


K 


2.2 


0.02 


-2.18 




L 


3.5 


0.01 


-2.51 




M 


4.9 


0.01 


-2.23 




N 


12.0 


0.0 


-3.58 





(1) Johnson et al. 
(3) Van Eck et al. 



(1966); (2) Thomas et al. (1976); 
(2000b); (4) Smith (2003) 



3. Simultaneous fit of observed visibilities and SED 

3.1. Photometric observations of it 1 Gruis 

Besides the ISO/SWS spectrum (see Fig. 1), ground-based 
photometry is also available for n 1 Gruis (see Table 2). The 
K magnitude has already been dereddened by Van Eck et al. 
(2000b). The other magnitudes were dereddened using the 
Fitzpatrick absorption curve for the M-band (Fitzpatrick, 
1999) and the values of Rieke & Lebofsky (1985) for the 
N-band. Van Eck's visual extinction value of A v =0.17 was 
adopted, so that the ratio of total to selective extinction 
at V, 3.1 on average, allowed us to correct the other mag- 
nitudes. At wavelengths longer than 12.5 /um, the inter- 
stellar extinction becomes negligible compared to the mea- 
surement precision of the magnitude. Table 2 lists the basic 
photometric properties of tt 1 Gruis. 

3.2. The stellar atmosphere of n 1 Gruis 

The realistic MARCS code (see Appendix A.l) was used 
to determine the synthetic intensity distribution of the 
S-star 7T 1 Gruis. Such a distribution is far more realistic 
than a blackbody distribution due to the large number of 
absorption lines in S star atmospheres. 
The 7T 1 Gruis MARCS model atmosphere is taken from a 
grid specifically designed for S stars (Plez et al. (2003) 
and in preparation), since their non-solar C/O ratio and 
s-process abundances impact their atmospheric structure. 
Synthetic J — K and V — K colors were computed for 
all models from this grid (Fig. 5). These reveal that the 
C/O ratio has a strong influence on the colors, especially 
for low effective temperatures (< 3200 K). There is, 
however, hardly a one-to-one relationship between the 



Table 3. Stellar atmospheric parameters of it 1 Gruis, all from 
Vanture & Wallerstein (2002) except C/O. 



Parameter 


Value 


T c ff 


3000 K 


log g 


0.0 


Wmicro (kmS _1 ) 


3.0 


[Fe/H] 


0.05±0.48 


[Zr/Ti] 


1.27±0.23 


C/O 


0.97 



photometric indices and the atmospheric parameters, 
and several models are found to match the (J — K) , 
(H — K)o and (V — K) indices of n 1 Gruis. All these 
models have T cff = 3100 K, log 3 = 0, and [Fc/H] = 0, 
but allow some range in [s/Fc] and C/O (namely C/O 
= 0.97 for [s/Fe] = 0.5 dex, and 0.90 < C/O < 0.92 for 
[s/Fe] = 1 dex). The situation is illustrated in Fig. 5, 
where n 1 Gruis falls within the dashed box. It must be 
noted that Vanture & Wallerstein (2002) have performed 
a spectroscopic abundance analysis of tt 1 Gruis, though 
not with dedicated MARCS model atmospheres, and 
found [Fe/H] = 0.05 ± 0.48 and [Zr/Ti] = 1.27 ± 0.23 
for T e fj = 3000 K. These are in good agreement with the 
parameters inferred from the MARCS models. Table 3 
summarizes the atmospheric parameters and the most 
important abundances adopted for n 1 Gruis. 



3.3. Analytical modeling 

Now that the MARCS intensity distribution of tt 1 Gruis is 
known, it can be used to model the spectro-interferometric 
measurements of the star. 
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Fig. 5. Dereddened (V-K,J-K)o color-color diagram for S stars 
(open triangles) from Van Eck & Jorissen (2000a) and Van 
Eck et al. (2000b), M giants (open squares), and a few C stars 
(crosses) both from Whitelock et al. (2000), compared to syn- 
thetic MARCS model colors (filled color symbols; Plez et al. 
2003) of different temperatures: 4000, 3800, 3600, 3400, 3200, 
3000 K. The black solid curve from lower left to upper right is 
the line joining all models with C/O ratio of 0.97. The ir 1 Gruis 
color range is enclosed in the dashed box. 

3.3.1. The analytical thin-dusty-shell model 

The general expression of the analytical thin-dusty-shell 
intensity distribution was derived by Cruzalebes & Sacuto 
(2006). These authors present an analytical solution to the 
radiative transfer equation in spherical geometry for an op- 
tically thin dusty shell. The model, briefly explained below, 
was compared to the well-established radiative transfer 
code DUSTY used to fit the spectro-interferometric data 
of 7T 1 Gruis (see Sect. 3.4.1). Using this thin-dust-shell 
model is justified by the low mass-loss rate of the star 
(4.6xl0~ 7 M Q yr -1 ; Jorissen & Knapp 1998) in view of the 
low-contrast, mid-infrared emission features. This led Gail 
(1990) to propose that the dust formation is insufficient 
for triggering a dust-driven outflow. The advantage of this 
model is the straightforward way a first set of physico- 
chemical parameters can be determined along with their 
errors. These output parameters are then inserted into the 
widely-used numerical radiative transfer code DUSTY in 
order to confirm and complete the analysis of this star (see 
Sect. 3.4). 

The expression of the intensity given by the analyt- 
ical thin-dusty-shell model is described in Sect. 2.2 in 
Cruzalebes & Sacuto (2006), where the central star is repre- 
sented by the corresponding synthetic MARCS S-spectrum 
(see Sect. 3.2). The global optical depth related to the in- 
frared continuum excess is given by 



Ao J V £ out 



where re\ n is the extinction coefficient at the inner 
boundary at the reference wavelength Ao- The extinction 
coefficient is parametrized as a decrease with wavelength 
as A _/5 . Values of (3 are usually taken between 0.2 and 
2.4 depending on the chemistry, nature (amorphous or 
crystalline), shape, and size of grains (Dupac et al., 2003). 
The factor p corresponds to the radial number density 
falling off as r~ p . The parameters £ m and e out are the 
angular radii of the shell inner and outer boundaries, 
respectively. Finally, Z out (e) is the total geometrical depth 
of the shell. 

The formalism described in Cruzalebes & Sacuto (2006) 
applies in the context of the infrared continuum excess 
caused by the dust shell. However, as we can see in Fig. 1, 
the spectrum of ir 1 Gruis exhibits a bump around 11 /<m 
caused by the emission from oxygen-rich dust. The opti- 
cal depth of this bump can be modeled with a Gaussian 
distribution (Groenewegen et al., 1994), 



bump 



A e 



-(^) 2 



(4) 



where A is the strength of the bump, A c the central wave- 
length of the dust feature, and AA related to the width 
of the bump. Therefore, the expression for the intensity 
must depend not only on the optical depth related to the 
infrared continuum excess (Eq. 3) but also on the optical 
depth related to the mid- infrared bump (Eq. 4). In the case 
of spherical symmetry, the visibility function V\ measured 
by an interferometer is linked to the angular distribution 
of the intensity of the source I\ (e) through the normalized 
Hankcl transform, 



h(Q) 



(5) 



where bij is the modulus of the projected baseline (on the 
plane of the sky) formed by the pair of apertures and 



h (f) = / 2neh (e) Jo (2nfe) de 



(3) 



(6) 



is the Hankel transform of I\ at the spatial frequency / = 
bij/X, where Jo is the zeroth-order Bessel function of the 
first kind. Note that /a(0) corresponds to the flux density 
of the source F\. 

The aim of this model is to give a first set of physico- 
chemical parameters by simultaneously fitting spectromet- 
ric and interferometric data. These parameters are then 
used as input parameters in the numerical code DUSTY 
to determine the final set of parameters characterizing the 
dusty circumstellar environment of 7T 1 Gruis. 

The independent parameters of the models can be deter- 
mined using the well-known iterative Levenberg-Marquardt 
algorithm minimizing the \ 2 function obtained by compar- 
ing modeled and measured visibilities, using the weight as- 
sociated with the N visibility measurements. These weights 
are determined by the inverse square of the visibility errors. 
The minimization process uses the analytical expressions 
for the first derivatives of the visibility function with respect 
to each independent parameter as extensively described in 
Cruzalebes & Sacuto (2006). 
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3.3.2. Spectro-interferometric results 

In this section, the ir 1 Gruis interferometric data are an- 
alyzed by means of the analytical model described in 
Sect. 3.3.1. The fit is done separately on two data sets, re- 
ferred to as DATA_SPECTRO.CAL and DATAJNT.CAL 
in the following. The former set corresponds to the ir 1 Gruis 
MIDI data calibrated with the 26.8 ± 1.3 mas diameter 
of [3 Gruis derived from the spectrometric method (see 
Appendix A.l). The latter set is calibrated with the 28.8 ± 
0.6 mas diameter of j3 Gruis derived from the interferomet- 
ric method (see Appendix A. 2). 

Some parameters of the analytical model were fixed at 
pre-determined values or only allowed to vary within a rea- 
sonable range: 

— The effective temperature (T c ff), which was found to 
be 3000 to 3100 K both from the work of Vanture & 
Wallcrstein (2002) and from the position of the star 
in the (J-K,V-K)o color-color diagram. Therefore both 
3000 and 3100 K synthetic MARCS models were utilized 
for the fit. 

— A dust temperature at the inner boundary (Tj n ) be- 
tween 600 and 800 K, in accordance with the low mass- 
loss rate of this star for which dust condensation will 
halt at an intermediate stage of the condensation-ladder 
(see Sect. 3.4.2). This is called the "freeze-out" mecha- 
nism (Heras & Hony, 2005; Blommaert et al., 2006). 

— The central wavelength of the dust feature (A c ), which 
is fixed to 11 /xm from the position of the mid-infrared 
bump excess (see Fig. 1). 

— Varying the density power-law coefficient (p) between 1 
and 2.5. 

— Varying the value of the spectral index (/?) between 0.2 
and 1.4, in agreement with values observed in circum- 
stellar environments (Weintraub et al., 1989; Beckwith 
et al., 1990; Knapp et al., 1993). 

The model provides in both cases good quality fits 
ix 2 /(N — p) = x^ od <0.6, where N is the number of mea- 
surements and p is the number of adjustable parameters); 
however, the parameter values are relatively different for 
the two data sets, especially for the central star diameter. 
Because this parameter is the most constraining for the 
stellar spectrum, one can easily distinguish which of 
the parameter sets provides the most reliable fit to the 
spectrum of the star by superimposing the analytical 
flux model generated with the 2 parameter sets onto the 
spectro-photometric data of the star. Figure 6 compares 
the 7T 1 Gruis fluxes predicted by the analytical model using 
the parameters for the data sets DATA.SPECTRO.CAL 
or DATAJNT.CAL, with the MIDI flux, the ISO/SWS 
spectrum, and the photometric data. It is very clear from 
this figure that the parameter set providing the best fit 
to the spectrum is the DA TA _SPECTR 0_ CA L one (see 
Table 4). 

The stellar diameter obtained from the fit of the analyt- 
ical model is </>=21.6 mas. This corresponds to R*=357 R© 
for the adopted Hipparcos distance of D=153 pc. In a study 
of the angular size of a sample of S and C stars, Van Belle 
(1997) found that, for non-Mira stars, S-type stars are 
on average 110 Rq larger than M-type stars. Since the 
relationship between radius and V — K color derived by 



Table 4. Parameters of tt 1 Gruis deduced from the fit of the ana- 
lytical visibility model to the data sets DATA.SPECTRO.CAL. 



Parameter Value 

Effective temperature (K) : T cft 3000 

Luminosity (L©) : L* 9300±2700 

Central star diameter (mas) : 2 21.6±0.3 

Shell inner radius (mas) : £i n 159±10 

Inner boundary temperature (K) : Ti n 776±25 

Dust bump wavelength (/im) : A c 11* 

Bump amplitude (xl0~ 4 ) : A 8.3±0.6 

Bump width (^m) : AA 2.8±0.3 

Bump optical depth (xl0~ 3 ) : ru Mm 1.5±0.2 

Density power law coefficient : p 1.7 

Spectral index : (3 1 

xLi 0-55 

* : fixed at pre-determined value. 



Van Belle (1999) applies to M stars, the corresponding ra- 
dius of 7T 1 Gruis is set to 359 R©, a value very close to the 
one obtained from the fit of the analytical model. In Table 4, 
the luminosity of the central star has been added. This one 
is deduced from the effective temperature of the central star 
(T cff =3000 K), its angular diameter (^=21.6±0.3 mas), 
and its distance (D=153±20 pc), 

L * = A7T {y D ) aT ° a ~ ( 9300± 270 °) L o> ( 7 ) 

where a is the Stcfan-Boltzmann constant. The correspond- 
ing apparent bolometric magnitude of n 1 Gruis is 0.8±0.6, 
in agreement with the value of 1.07 found by Van Eck et 
al. (2000b). 

Another significant parameter value is the density power- 
law coefficient. With an r~ 17 density distribution, this sug- 
gests a decrease in mass-loss rate with time, in line with 
the finding of Chiu ct al. (2006) that the mass-loss rate 
of tt 1 Gruis has probably decreased over the past 90 years. 
This set of parameters will thus be used as input for the 
DUSTY model described in Sect. 3.4.1. 

3.4, Numerical modeling 
3.4.1. The DUSTY code 

The DUSTY code (Ivezic et al., 1999) was used to generate 
synthetic spectrally-dispersed (8-13 /im) flux and visibilities 
to confirm and complete the parameter set found with the 
analytical model (see Table 4). DUSTY is a public-domain 
simulation code that solves the problem of radiation trans- 
port in a circumstellar dusty environment by integrating 
the radiative transfer equation in plane-parallel or spheri- 
cal geometries (Ivezic & Elitzur, 1997; Ivezic et al., 1999). 

The DUSTY model integrates the MARCS synthetic 
stellar spectrum of the central star (see Sect. 3.2) and con- 
siders as input parameters a part of those deduced from the 
fit of the analytical model: T e ff, L*, Ti n , th^, and p (see 
Table 4). DUSTY also requires the grain size distribution 
as input so we chose a standard MRN grain size distribu- 
tion as described by Mathis et al. (1977): n(a) cx a~ 35 with 
minimum and maximum grain sizes of 0.005 /im and 0.25 
/im, respectively. 

Therefore, the analytical model allows one to reduce 
the number of DUSTY input parameters from 8 to 2 
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Fig. 6. Comparison of the n 1 Gruis fluxes predicted by the 
analytical model using the parameters for the data sets 
DATA.SPECTRO.CAL (solid line) or DATAJNT.CAL (dashed 
line), with the MIDI flux (error bars), the ISO/SWS spectrum 
(blue line), and the photometric data (open triangles). 



(namely the chemical composition of the dust and the 
thickness of the shell). Indeed, the chemical composition 
of the dust is related to the value of the spectral index 
p. However, this value encapsulates many characteristics 
related not only to the chemistry of the grains but also to 
their nature (amorphous or crystalline), their shapes, and 
their sizes. 

The thickness of the shell was found to be the least 
accurate output parameter, with a lower limit of 60 stellar 
radii, corresponding to the beam size of an AT. In fact, 
varying this value from 60 to 1000 times the stellar radius 
does not have a significant impact on the model. This is 
related to the low mass-loss rate of it 1 Gruis (4.6xl0~ 7 
MQyr -1 ) leading to a very low dust density in the outer 
layers. Since the N band is not well-suited to determining 
the outer radius of the dust shell, this one could extend to 
the inner radius of the flared disk located at 125 stellar 
radii (Chiu et al., 2006). 



3.4.2. Dust chemical composition 

The last free parameter to determine for the DUSTY mod- 
eling is the chemical composition of the dust shell. However, 
this parameter can be constrained from the shape of the 
mid-infrared spectrum. 

The classification by Little-Marenin & Little (1988, 
1990) (hereafter referred to as the LML classification) gives 
a sequence from the classic, narrow-emission feature at 
10 /zm (classified as silicate) to features with progressively 
stronger contributions at 11 /im (silicate+ and silicatcH — h, 
depending on its strength) to broad, low-contrast emission 
peaking longward of 11 /im (Broad). In this classification, 
Little-Marenin & Little (1988) investigate the LRS (IRAS 
Low-Resolution Spectrometer) spectra of a large number of 
MS and S stars and suggest that the unusual dust chem- 
istry (caused by a C/O ratio near unity) might produce 



unusual dust excess. Little-Marenin & Little (1988) create 
a separate class of silicate emission, the S class, which they 
associate with S stars and describe as having a peak emis- 
sion around 10.8 /im, which is the case for tt 1 Gruis (see 
Fig. 1). However, later analysis by Chen & Kwok (1993) 
and Sloan & Price (1995) concluded that the LRS spectra 
of S stars are not significantly different from those of M 
giants. Sloan & Price (1998) noticed that MS and S stars, 
having a C/O ratio near unity, produce circumstellar shells 
dominated by alumina grains. They suggest that, when the 
C/O ratio approaches unity, the formation of CO leaves 
little oxygen available for grain formation. As gaseous ma- 
terial moves away from the central star and cools, alumina 
grains will form before silicate grains because of the higher 
condensation temperature of alumina. Sloan & Price (1998) 
also suggest that stars with C/O ratios near unity may ex- 
haust the supply of available oxygen with the formation of 
CO gas and alumina grains, preventing silicates from form- 
ing. 

Figure 1 shows the broad dust emission feature from 9 to 
13 /im peaking at 10.5 /im, which originates from the opti- 
cally and geometrically thin shell of tt 1 Gruis. This spectral 
shape is well-known and corresponds to one phase of the 
silicate dust sequence scenario (hereafter referred to as the 
SDS scenario) introduced by Stencel et al. (1990). In this 
scenario, grains form first from condensing alumina ma- 
terial. Then, as the grains evolve, silicate material begins 
to dominate the emission. Fresh silicate material would be 
crystalline and exhibit rather narrow emission components 
at 10 and 11 /im, while grains in more extended shells would 
be amorphous and exhibit the classic silicate feature. The S 
star 7T 1 Gruis belongs to the bottom end of the SDS, classi- 
fied as SE2 by Sloan & Price (1998). This lower silicate dust 
sequence is populated by spectra with broad, low-contrast 
emission features, originating from optically thin shells of 
alumina dust. 

In this context, we first tried to reproduce the mid-infrared 
spectral shape emission with a pure alumina shell using 
porous amorphous AI2O3 of Begemann et al. (1997) 2 re- 
trieved from the Database of Optical Constants of the 
Laboratory Astrophysics Group of the AIU Jena 3 . This 
composition failed to fit the spectrometric data. Spectra of 
SE2 class stars do not exhibit any emission peak at 10 /im. 
However, the ISO/SWS and MIDI spectra of tt 1 Gruis show 
that the emission feature begins to peak at 10 /im, ending 
beyond 12 /im (see Fig. 1). Therefore, in spite of the SDS 
classification of tt 1 Gruis as SE2, it seems that this star 
is close to the structured silicate emission class (SE3-SE6). 
This class is intermediate between the silicate and the broad 
class stars implying that the model contains a mixture of 
silicate and alumina grains, as discussed by Lorenz-Martins 
& Pompeia (2000). Consequently, we decided to work with 
mixtures of aluminum oxide and warm silicate (Ossenkopf 
et al., 1992), ranging from 100% alumina grains to 100% 
silicate grains in 10% increments. We find a good fit (see 
Fig. 7) with a ratio of 30% alumina and 70% silicate grains 
(see Table 5). 



2 Note that the data of Begemann et al. (1997) only cover 
emission between 7.8 and 500 /zm. DUSTY assumes a constant 
refractive index in the unspecified shorter wavelength range, 
with a value equal to the corresponding first point of the tabu- 
lation. 

3 http : //www. astro .uni-jena. de/Laboratory/Dat abase/ 
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Table 5. Parameters for tt 1 Gruis deduced from the fit of the 
DUSTY model to the interferometric and spectro-photometric 



data. 


Parameter 


Value 


Effective temperature (K) 


3000+ 


Luminosity (Lq) 


8700+ 


Central star diameter (mas) 


21.0 


Shell inner radius (mas) 


149.5 


Inner boundary temperature (K) 


700+ 


Grain chemical composition 


70% WS + 30% Al. 1 


Density power law coefficient 


1.7+ 


Grain size distribution 


MRN* 


Geometrical thickness (e*) 


60 2 


ff fim optical depth (xfCT 3 ) 


8.3+ 


* : fixed at pre-determined value. 



+ : derived from the analytical model (see Table 4). 
1 WS stands for warm silicates and Al. for alumina. 
2 This value is not well defined (see text). 



The 11 ^m global optical depth found with the analyt- 
ical method (rn Aim =(1.5±0.2) x 1CT 3 ; Table 4) appears to 
be somewhat too small for DUSTY, which requires about 5 
times more optical depth to yield a good fit of the IR spec- 
trum. This may be due to the analytical model assuming 
that the 11 /xm feature may be approximated by a Gaussian 
function, whereas Fig. 1 shows that this excess is clearly 
asymmetric. 

The temperature at the inner boundary of the 
shell (700 K; Table 5) is below the aluminum oxide 
(T^° 3 «1500-1800 K) and silicate (T^ nd wl000-1500 K) 
condensation temperatures. Such low temperatures in the 
shells of evolved stars were first suggested by Rowan- 
Robinson & Harris (1982) and later by Onaka et al. 
(1989) and Simpson (1991), based on studies of the 
IRAS/LRS data. Onaka et al. (1989) argue that these 
low temperatures cannot be understood in the framework 
of homogeneous nucleation theories and propose, as a 
possible explanation, that mantle growth on pre-existing 
aluminum oxide particles is the major process of sili- 
cate formation. Extending the temperature at the inner 
boundary of the shell to alumina condensation temper- 
ature (^1500 K) can also explain both the spectrum 
and visibility of the star. However, assuming a mass-loss 
rate value of 4.6xl0~ 7 M Q yr _1 for it 1 Gruis (Jorissen & 
Knapp, 1998) and a dust-particle specific density of 3 
g cm~ 3 , the value of the gas-to-dust ratio of the model 
is 60. We obtain a value of 165 for the inner boundary 
temperature of 700 K, which agrees with the gas-to-dust 
ratio of 158±13 obtained by Knapp (1985) for oxygen stars. 

Figure 7 shows the best interferometric and spectro- 
photometric fit of the DUSTY model corresponding to the 
parameter set of Table 5. The model appears to match the 
observed spectra well and the MIDI visibilities fairly well. 
However, while the spectral fit of the DUSTY model is good 
overall but especially so in the mid-infrared band (8-13 /im), 
the model visibilities are over-evaluated for the two ^15 
and 30 m projected baselines in the 8-9 /ttm spectral band 
(see Fig. 7). This reveals that the model is not extended 
enough for these baselines at these wavelengths. This can 
be explained by the very low optical depth in that spec- 
tral range (f (8-9 /zm)=2.5x 10~ 3 ) leading to a domination 
of the central star flux over the flux of the dust shell. In 



the 8-9 /im spectral range, the model therefore has values 
very close to the central star's uniform-disk intensity dis- 
tribution, leading to an increase in the visibility. Moreover, 
one can see that the 8-9 /im visibilities decrease with the 
increase in the baseline lengths, which is understandable 
in the context of an equivalent uniform-disk visibility func- 
tion. 

A way to increase the spatial extension of the model be- 
tween 8 and 9 /zm is to add a molecular shell between the 
photosphere of the star and the inner boundary of the dust 
shell. In the next section, we investigate the possibility of 
adding such a molecular shell in the model. 

3.5. Addition of a warm H20+SiO molecular shell 

Recent work reveals that diameters of the Mira and non- 
Mira stars appear systematically larger than expected in 
the near-infrared (Mennesson et al., 2002; Ohnaka, 2004b; 
Perrin et al., 2004; Weiner, 2004) and in the mid- infrared 
(Weiner et al., 2000, 2003a,b; Weiner, 2004; Ohnaka, 2004b; 
Ohnaka et al., 2005). Such an increase cannot be attributed 
only to dust shell features, but also to the possible existence 
of extended gaseous shells close to the stellar photospheres. 
The analyses of the spectra obtained with ISO/SWS 
have revealed that the spectra of Mira variables are well- 
reproduced by water vapor shells extending to 2-4 stellar 
radii with gas temperatures of 1000-2000 K and column 
densities of 10 20 -10 22 cm~ 2 (Tsuji et al., 1997; Yamamura 
et al., 1999; Matsuura et al., 2002). Since tt 1 Gruis is a non- 
Mira star, we should expect different parameter values than 
those for Mira variables. 

The single dust shell modeling with DUSTY (see Fig. 7) 
suggests that we need to have a contribution from a molec- 
ular shell only in the 8-9 fxm spectral band. Such a molec- 
ular shell, also called "MOLsphere" (Tsuji et al., 2000), is 
favored by the low surface gravity (log g=0; see Table 3) and 
high luminosity (^9000 L@) of the star, allowing the levita- 
tion of the upper layers. The temperatures in this levitated 
matter can be low enough for molecules to form (Verhoelst 
et al., 2006). 

Because this molecular shell has a larger extent than 
the central star, the former is resolved with smaller 
baselines so as to decrease the visibility amplitudes in 
that spectral range. At the same time, the molecular 
shell must make a negligible contribution to the total 
flux in the 9 to 13 /im spectral range relative to the 
dust shell, to recover the good fit that we found with 
the single dust shell model at these wavelengths (see Fig. 7) . 

Water vapor (H2O) and silicon monoxide (SiO) are the 
main molecules contributing to absorption and causing the 
apparent mid-infrared diameter to be larger than the con- 
tinuum diameter (Ohnaka et al., 2005; Verhoelst et al., 
2006; Perrin et al., 2007). The H2O lines are present across 
the entire N band range, while most SiO lines are located 
between 8 and 10 ^m. 

The model of molecular shell that we chose is the same 
as the one used by Ohnaka et al. (2005) and consists of 
H2O and SiO gas with constant temperature and density, 
extending from the central star radius up to the external 
boundary e mo i. The input parameters of the model are 
the temperature of the molecular shell (T mo i), the column 
densities of H 2 and SiO (N^ and N^p, respectively) in 
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Fig. 7. Top-left: best fit of the dust shell model spectrum (red line) on the MIDI flux (error bars), the ISO/SWS spectrometric 
data (blue line), and the photometric data (open triangles) of tt 1 Gruis. The central star contribution (dotted line) and the dust 
shell contribution (dashed line) are added. Top-right: close-up view of the best fit of the dust feature. Bottom: the corresponding 
model visibility (solid line) superimposed on the MIDI visibilities (error bars) for the seven projected baselines. 



the radial direction, and the geometrical extension of the 
shell (e mo i). 

In this model, we first calculate the line opacity due to 
H2O and SiO with a wavenumber interval of 0.04 cm -1 
over a wavelength range between 8 and 13 /im. We adopt 
a Gaussian line profile with an FWHM of 5 kms -1 , which 
represents the addition of a gas thermal velocity of a 
few kms" 1 (V T =1.5(T mol /300 K) 05 kms" 1 ) and the 
micro-turbulent velocity in the atmosphere of tt 1 Gruis 
(3 kms -1 ; see Table 3) and assumes that the molecular gas 
is in local thermodynamical equilibrium. The list of H2O 
lines was taken from the HITEMP database (Rothman, 
1997) while that of the fundamental bands of SiO was 
generated from the Dunham coefficients given by Lovas 
et al. (1981) and the dipole moment matrix elements 
derived by Tipping & Chackerian (1981). 

The resulting intensity distribution of the object is com- 
posed of the central star with a MARCS model atmosphere, 
the molecular shell, and the dust shell. The analytical ex- 
pression of the model is given by the expressions below: 

im = /r rcs e-— e"^ n (JL\ (8) 
IT°\e) = B x (T mol ) e"— [l - e —">] II (9) 



Table 6. Best combination of input values for the determination 
of the molecular shell. 



Component T mo i (K) e mo i (e*) 


N co i (cm- 2 ) 


H2O 1000 4.4 


5xl0 18 


SiO 


2xl0 18 



I^ st (e) = B x (T dust ) [1 - ] n f_j (10) 

where r mo i and Tdust are the optical depths of the molecu- 
lar and dust shells, respectively, and IT is the uniform disk 
distribution (II(x) = 1 if —1/2 < x < 1/2, zero otherwise). 
The corresponding flux and visibility, obtained by the 
Hankel transform of the intensity profile, are then spec- 
trally convolved to match the low spectral resolution (30) 
of the MIDI prism. Using the above model and without 
changing the parameter values found for the thin dust shell 
(see Table 5), we search for the best combination of in- 
put parameters (T mo i, e mol , N co 2 j , and N^°) to improve 
the spectro-interferometric fitting in the 8-9 /im spectral 
range. The best-fitting values are presented in Table 6. 

These values correspond to an extended, optically thin 
(T m oi<0.1) molecular shell. Such a diluted molecular shell 
can be connected to the nature of tt 1 Gruis. First, because 
the C/O ratio is close to unity for S stars (C/O=0.97 for 
7T 1 Gruis; see Fig. 5), the last oxygen-bearing molecules dis- 
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appear (Ferrarotti & Gail, 2002). Second, due to its classi- 
fication as a non-Mira semi-regular variable, for which the 
variability amplitudes are much smaller than for Mira vari- 
ables, a less dense molecular shell is indeed expected. From 
these arguments, it seems reasonable that this star exhibits 
weak spectral features from the molecular shell. 

Figure 8 shows our best spectrometric and interferomct- 
ric fit to the optically thin H 2 0+SiO molecular shell and 
optically thin silicate+alumina dust shell model using the 
parameter sets of Tables 5 and 6. 

This figure reveals that the spectrometric, as well as 
the visibility measurements are now well-reproduced over 
all the 8-13 fim spectral band for all projected baselines. 
In particular, the model visibilities have decreased in 
the 8-9 (im spectral band compared to the single dust 
shell model (see Fig. 7). This is due to the correlated 
flux contribution of the molecular shell leading to a more 
extended object at these wavelengths. With a size of 
4.4 stellar radii, the molecular shell is almost completely 
resolved with baselines larger than 30 meters. For this 
reason, its contribution to the correlated flux is greater 
for the 15- and 30-meter projected baselines than for the 
60-meter one. The optical depth of the molecular shell 
outside the 8 to 9 /j,m range is ten times less than in this 
band. Therefore, the correlated flux contribution of the 
molecular shell is almost negligible at these wavelengths, 
preserving the good fit from the single dust shell model. 

From this last model, we can now infer the spatial cir- 
cumstellar structure of it 1 Gruis. It consists of an optically 
thin warm H^O+SiO molecular shell extending from the 
stellar photosphere out to 4.4 stellar radii, then an opti- 
cally thin alumina-silicate dust shell beginning at about 14 
stellar radii, and finally a slowly expanding flared disk with 
a central cavity of radius about 125 stellar radii, determined 
by Chiu et al. (2006) in the millimetric spectral range. 

4. Asymmetries of the circumstellar material and 
binarity 

The bipolar outflow perpendicular to the flared disk was 
discovered by Sahai (1992) and subsequently studied by 
Knapp et al. (1999a) and Chiu et al. (2006). A similar 
structure combining disk and jet was observed by Knapp 
et al. (1997) in another very evolved AGB star, VHya, 
whose binary nature was subsequently discovered by 
Knapp et al. (1999b). As argued by Huggins (2007), the 
binarity seems to be a necessary condition creating a 
disk/jet structure in evolved mass-losing stars. Although 
7T 1 Gruis has a GOV companion located 2.71" away, this 
companion falls within the expanding disk (Chiu et al., 
2006), a strange situation that has not been encountered in, 
for example, binary post-AGB stars, where the companion 
falls well within the cavity interior to the disk (Waclkcns 
et al., 1996; Dominik et al., 2003). Therefore, Chiu et al. 
(2006) wonder whether there 'might be a much closer 
companion that has so far escaped detection and currently 
is the source of the high- velocity outflow'. 
Makarov & Kaplan (2005) and Frankowski et al. (2007) 
find possible evidence for such a close companion from 
the discrepancy between the Hipparcos and Tycho-2 
proper motions. The Hipparcos proper motion relies on 
measurements only spanning 3 years, whereas the Tycho-2 



one uses positional data spanning more than a century. 
Any discrepancy between these short-term and long-term 
proper motions is thus indicative of an unrecognized 
orbital motion in addition to the true proper motion. The 
discrepancy is very significant in the case of ir 1 Gruis. 
According to the analysis of Frankowski et al. (2007), these 
so-called 'A/i binaries' correspond to systems with orbital 
periods in the range 1500-10 4 days, much shorter than the 
6000 years period inferred by Knapp et al. (1999a) for the 
GOV companion 2.7" away. 

We now estimate the likely range for the orbital separa- 
tion of the putative close companion and compare it to the 
size of the various shells and to the stellar radius. The pe- 
riod range for A/i binaries with known orbital elements are 
from 1500 to 10000 days, according to Fig. 4 of Frankowski 
et al. (2007). Adopting a total mass of 2 M© for the sys- 
tem, the corresponding range in relative semi-major axes 
are 3.2 to 11.4 AU, or 20.8 to 74.4 mas (for the parallax of 
6.5 mas). Adopting 3 M Q instead yields 3.7 to 13.1 AU or 
24 to 85 mas. These values must be compared to the 

— stellar radius = £ sta r = 10.5 mas 

— H20+SiO molecular radius = 4.4 e star = 46.2 mas 

— inner dust shell radius = 14.2 e s t a r = 149.5 mas. 

It is seen that only the longest periods and highest total 
masses are compatible with the shell sizes. One can wonder 
what in this case would be the impact of a companion 
orbiting at about 80 mas from the central star on the 
molecular and dust shell. As the extended dust shell 
is over-resolved, only the molecular shell could provide 
some hints of the presence of a perturbation. However, 
the mass-loss rate is so low that a high dynamic range 
is required to detect the asymmetries from the faint 
contribution from the molecular and dusty material. Such 
a dynamic range cannot correctly be achieved with the 
current MIDI capabilities. Nevertheless, the constraints 
provided by the differential phases that do not exhibit any 
obvious departure from spherical symmetry (see Fig. 3) 
have to be considered, since it is more sensitive to any 
departure from centro-symmetry. 

The A/i method offers no way to infer the position angle 
of the companion. However, if it is powering the jet, its 
orbital plane is likely to be coincident with that of the 
flared disk. The MIDI intcrferometric uv coverage (PA=41 
to 102°) is oriented in the direction of the flared disk 
(and thus along the orbital plane), as deduced from Fig. 4 
in Chiu et al. (2006). It is then difficult to envisage the 
presence of a close companion without any detection of 
asymmetries. This could be the case if the dust shell were 
not heated by the companion, which must therefore be on 
the lower main sequence (cool and faint). 

Beyond there being no detection of asymmetries as 
shown by the zero-differential phase (see Fig. 3), there 
is a more fundamental information inferred from the 
intcrferometric observations in Fig. 4. There is no obvious 
indication that the star alone could drive the jet and 
shape its surrounding. Notwithstanding this fundamental 
constrain, we are left with two possibilities, each of them 
invoking the influence of a companion: 

(i) There is indeed a close companion (cool and faint) 
to power the jet, but then why are there no departures 
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Fig. 8. Top-left: best fit of the proposed HbO+SiO molecular shell and silicate+ahimina dust shell model spectrum (red line) on 
the MIDI flux (error bars) and the ISO/SWS spectrometric data (blue line) of tt 1 Gruis. The central star contribution (dotted line), 
the molecular shell contribution (dot dashed line), and the dust shell contribution (dashed line) are added. The other figures show 
the corresponding model visibility (solid line) superimposed on the MIDI visibilities (error bars) for the seven projected baselines. 



from spherical geometry close to the star, related to the 
presence of an accretion disk? 

(ii) The close A/i companion is spurious, but then 
the jet must be powered by the far-away companion, and 
why is this outer companion orbiting the disk without 
creating a gap in it, as would be evident in the millimctric 
observations? 

Regarding item (i), one explanation could come from 
the close, optically thin environment of tt 1 Gruis leading 
to similar SEDs for accretion disks and spherical dust 
shells, as discussed by Knapp et al. (1993) and van der 
Veen et al. (1994). If we add to this argument the fact 
that the uv coverage is preferentially oriented in the 
direction of the disk having an opening angle larger than 
the interferometric angular coverage (^60°), it is then easy 
to confuse this structure with a spherical shell. 

Regarding item (ii), it must be mentioned that the A/i 
method might face problems when dealing with highly 
evolved stars and might then yield spurious results. To 
substantiate this, notice that Table 8 of Frankowski et al. 
(2007) lists several other candidate A/x binaries among 
long-period variables and suggests that these detections 
are similarly spurious, although the reason thereof remains 
unidentified. 

Still, even in the absence of a close companion, the 
Huggins (2007) analysis of the jet /torus properties in 
several AGB stars (among which tt 1 Gruis) and proto- 
planetary nebulae requires them to be binaries. The 2.7" 
companion to tt 1 Gruis must then be the one responsible 
for driving the jet/torus structure. In that case however, 
the conclusion of the CO data analysis by Chiu et al. 
(2006), which locates the companion within the flared 
disk, is quite surprising (when compared to situations 
prevailing in post- AGB systems for instance). Such a 
situation cannot be stable in the long term since the 
companion should carve a gap in the disk. If correct, 
it indicates that the disk is quite young, as confirmed 



from its kinematics. Taking into account a mean value of 
11 kms -1 for the wind expansion velocity found by Chiu 
et al. (2006) and a distance of 153 parsec (Perryman et 
al., 1997), we can estimate the typical transit time of the 
dust shell to be about 10 years, while the flared disk age 
is about 350 years (Huggins, 2007) with a strong decrease 
in the mass-loss rate about 90 years ago (Chiu et al., 2006). 



5. Summary and conclusions 

The goal of this work was to perform a multi-wavelength 
and multi-spatial resolution study of the close environment 
of 7T 1 Gruis with the N band spectral range of MIDI. Our 
new high spatial-resolution observations of the environ- 
ment of 7T 1 Gruis, consisting of MIDI/VLTI spectra and 
visibilities, combined with the existing millimeter data, 
led to a more comprehensive picture of its circumstellar 
environment on different spatial and spectral scales. 

A preliminary study of MIDI data shows that the 
mas-scale observations do not reveal any strong departure 
from the sphericity of the circumstellar environment, 
whereas arcsec-scale millimetric observations by Knapp 
et al. (1999a) and Chiu et al. (2006) show strong evidence 
of a slow and dense outflow (v~15 kms -1 ) and of a fast 
bipolar ejection (v>90 kms -1 ). Even though the con- 
straints provided by these interferometric observations are 
limited in terms of uv coverage and dynamics, this suggests 
that the companion located at 2.7" (400 AU at 153 pc, 
i.e. far outside the MIDI interferometric Field-Of-View of 
about 175 AU) may have some influence on the shape of 
the large-scale CO environment and may possibly launch 
the fast outflow that has been detected (cf. the models of 
Mastrodemos & Morris 1999 and the discussion in Chiu 
et al. 2006). 

The more detailed study starts with determining a 
synthetic atmosphere model for tt 1 Gruis to obtain approx- 
imate physical characteristics of the central star. The dusty 
circumstellar environment was initially analyzed with the 
optically thin analytical radiative transfer model developed 
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by Cruzalebes & Sacuto (2006). This model allowed us to 
determine some parameter values related to the central 
star and its dust shell, which we used to generate a more 
refined model with the numerical radiative transfer code 
DUSTY. This code helped us to determine the chemical 
composition of the close circumstellar environment, which 
was one of the main aims of this study. We show, from 
the literature, that this star belongs to one phase of the 
SDS scenario (Sloan & Price, 1998), in which grains first 
form from condensing alumina material. We confirmed this 
tendency using spectro-interferometric fitting and derived 
a ratio of alumina to silicate dust of 0.4, revealing that 
7T 1 Gruis is situated in an intermediate class named the 
structured silicate emission class. 

We then included an optically thin H 2 0+SiO molecular 
shell extending from the photosphere of the star up to 
4.4 stellar radii with a typical temperature of 1000 K 
to account for the fact that the single dust shell model 
is not extended enough in the 8-9 /zm spectral range. 
The formation mechanism of molecular shells in non-Mira 
semi-regular variables like tt 1 Gruis, in which the variability 
amplitudes arc much less than in Mira stars, can lead 
to less dense MOLspheres. Moreover, as S stars have 
a C/O ratio close to unity, neither oxygen-bearing nor 
carbon-bearing molecules are abundant in the atmosphere, 
with most of carbon and oxygen atoms locked up in CO 
molecules. 

Complex structures revealed by millimetric observations 
suggest that tt 1 Gruis may be at the end of its evolution on 
the AGB (Knapp et al., 1999a). This evolutionary stage 
is confirmed by the position of the star on the HR dia- 
gram of S stars (see Fig.7 of Van Eck et al. 1998). The 
large and diluted optically thin close environment around 
evolved stars could thus be part of the natural transition 
from the AGB to post- AGB and planetary nebulae stages. 
If this is true, what are the mechanisms responsible for 
such physical processes? Answers could come from long- 
term spectro-interferometric observations of the star and 
better Fourier uv coverage, which could be provided by the 
second generation VLTI instruments. This should be car- 
ried out in close synergy with ALMA array observations, 
which should provide images in the millimeter domain at 
similar angular resolution. 
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Appendix A: Angular diameter of the calibrator 
P Gruis 

The calibrator star (3 Gruis is a cool giant of spectral type 
M5111 without any detected circumstellar matter (Sloan & 
Price, 1998). With an estimated angular diameter around 
27 mas (Bedding et al., 1994; Heras et al., 2002), this star 
clearly cannot be considered as an unresolved object for the 
baselines used in this study. 

To determine the theoretical visibility of (3 Gruis, one 
can consider this star as having a uniform disk intensity. 
Such a hypothesis is valid in the mid-IR because of the 
vanishing limb-darkening coefficient at these wavelengths. 
At 10 /im, the ratio between the limb-darkened angular di- 
ameter and the uniform disk angular diameter is very close 
to one (Davis et al., 2000). The visibility of the uniform 
disk is 



'-if 



Jl (Wcal^ 



7T$cal : 



(A.l) 



where $ ca i is the angular diameter of the calibrator (3 Gruis, 
B p the projected baseline, A the observing wavelength, 
and Ji the first-order Bessel function of the first kind. 



Previous spectrometric (Heras et al., 2002) and interfer- 
ometric (Bedding et al., 1994) studies of this star allowed 
angular diameters of 28.1 mas and 27.0 mas, respectively, 
to be derived, although these evaluations bear an uncer- 
tainty of ±3 mas. This corresponds to an error bar on the 
equivalent uniform disk visibilities reaching 20% at high 
spatial frequencies. Moreover, this error comes on top of 
the transfer function error bar of 10% (see Sect. 2). Using 
the observed angular diameters of [3 Gruis would thus carry 
an uncomfortably large error bar of about 30%. We there- 
fore chose instead to predict the diameter of (3 Gruis by 
other means to reduce this error. Two methods are consid- 
ered in the following, one with a dedicated SED and the 
other a visibility model of the star. 



A.l. Spectrometric estimation 

A fit of a spectral energy distribution (SED) model to the 
spectrometric measurements of the calibrator allows the an- 
gular diameter of the calibrator (3 Gruis to be determined: 

2 

jcal = F cal ; (A. 2) 

where $ ca i is the angular diameter of the calibrator, 7^ al 
the intensity of the calibrator from a MARCS model, and 
F^ 1 the ISO/SWS spectrometric data flux of the calibrator 
[3 Gruis. 

We used the MARCS code (Gustafsson et al., 1975; Plez 
ct al., 1992, 1993), which generates stellar atmospheres of 
cool evolved stars, solving radiative transfer in a hydrostatic 
LTE environment using spherical geometry. To deduce the 
closest MARCS model for the calibrator (3 Gruis, we must 
know the stellar atmospheric parameters. These parameters 
have been derived by Judge (1986) using high-resolution 
spectra from the International Ultraviolet Explorer (IUE) 
satellite: 




- T cff =3400±100K 
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- log g=0.4 

- [Fe/H]=Q.O. 

The last stellar atmospheric parameter is the C/O 
ratio. The carbon abundance (log e( 12 C)=8.4) is scaled to 
the hydrogen (loge(H)=12) and oxygen (log e( 16 0)=8.73) 
average abundances for the sample of normal M giants 
derived by Smith & Lambert (1985). These lead to a C/O 
ratio of 0.47 for (3 Gruis. 

The U/B/V/R/I/K/L/M/[ll]/[22]/[60] photometric 
data are taken from Johnson et al. (1966); Feinstein (1966); 
Mendoza (1969); Engelke (1992). Dereddening is not nec- 
essary, since the galactic latitude of the star (-58°) leads 
to zero visual extinction (Chen et al., 1998). Errors on the 
angular diameter of the star can be estimated from the 
phases of the 2.36-45.39 /Ltm spectrometric data from the 
ISO/SWS (Infrared Space Observatory, Short- Wavelength 
Spectrometer) (Sloan et al., 2003). The spectra of this star 
were obtained on 1996 Apr 9 and 1997 May 7. Assuming 
a period of 37 days for the SRb variable with a V am- 
plitude of 0.34 mag, and using the ephemeris of Otero & 
Moon (2006), we find that the ISO observations are taken at 
phases 0.61 (i.e. shortly after the photometric maximum) 
and 0.99, at the photometric minimum. Thus, the MIDI 
phases of this star are approximately intermediate between 
the two ISO/SWS phases, except for the 2006 June 19 ob- 
servation corresponding to phase 0.01 (see Table 1). 

Figure A.l presents the best fits of the MARCS model 
spectrum on each of the two ISO/SWS spectrometric ob- 
servations (1996 Apr 9: Phase 0.61 and 1997 May 7: Phase 
0.99) yielding, from Eq. A. 2, the angular diameter of the 
calibrator. Both fits give an angular diameter of 28.0 mas 
for Phase 0.61 and 25.5 mas for Phase 0.99, leading to a 
mean angular diameter of 26.8 mas with a mean standard 
deviation of 1.3 mas. This error is about 2 times smaller 
than the one found in the literature (Heras et al., 2002) for 
(3 Gruis, thus reducing the corresponding error bar on the 
uniform disk visibility (from 0.6% at low spatial frequencies 
to 10% at high spatial frequencies). 
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Fig. A.l. Top: best fit of the MARCS model spectrum (red 
line) on the calibrator f3 Gruis ISO/SWS data (blue line) ob- 
tained on 1996 Apr 9 (Phase 0.61). Bottom: same as the 
previous one for 1997 May 7 (Phase 0.99). The photometric 
(U/B/V/R/I/K/L/M/[ll]/[22]/[60]) data of /3 Gruis (open tri- 
angles) are superimposed on each figure. 



A. 2. Interferometric estimation 

Another method for deriving the angular diameter of the 
calibrator (3 Gruis consists of using observations of several 
calibrators observed during the same nights but at a differ- 
ent sky position. The journal of observations for these stars 
(see Table A.l) shows their spectral type and their angular 
diameter with their errors calculated from the ESO VLTI 
Calibrators Program (Richichi & Percheron, 2005), with 
the position, UT time, and the relative projected baselines 
of the calibrated sources compared to the ones for (3 Gruis 
of each observation. 

This information illustrates the two main difficulties of 
these data. First, the time interval between the observation 
of (3 Gruis and its calibrator is rather long (from lh30 to 
4h30) in comparison to the mean intermediate time ded- 
icated to the observation of an interferometric calibrator 
star. The second difficulty concerns the difference in the sky 
position of these calibrators compared to that of (3 Gruis. 
The differences in position are rather large (from 1 to 4h30 
for the right ascension and from 10 to 55° for the decli- 
nation), which may also lead to markedly different atmo- 
spheric conditions between the source and calibration ob- 



servations. Consequently, the same difficulties that affect 
the determination of the calibrated visibility of it 1 Gruis 
affect the use of unresolved stars located at some distance 
on the sky for calibration. However, the atmospheric con- 
dition differences introduce only an amplitude bias on the 
calibrated visibility of the star (see Sect. 2.1), 

Vsou = (A.3) 

where a is a factor related to the atmospheric condition 
differences between the source and the calibrator observa- 
tions. In the case of a simple uniform-disk intensity dis- 
tribution like (3 Gruis, the equivalent angular diameter is 
independent of the a factor. This is not the case for more 
complex models that describe a source like ir 1 Gruis, which 
is surrounded by a circumstellar envelope. Such models use 
parameters related to the a factor, which can lead to incor- 
rect interpretation of the visibility measurements. 

Figure A. 2 shows the fits of the equivalent uniform disks 
to the calibrated spectrally-dispersed visibilities of (3 Gruis 
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Table A.l. Journal of observations of the calibrators observed the same nights as f3 Gruis. 



Star Spectral type Right ascension ; Declination UT date & Time duD[mas] Basc[m] PA[dcg] 



/3Gru 
ePeg 
/3Gru 
r?Sgr 
/3Gru 
t 4 Eri 



M5III 
K2Ib 
M5III 
M3.5III 
M5III 
M3/M4III 



22 42 40 

21 44 11 

22 42 40 
18 17 38 
22 42 40 
03 19 31 



-46 53 04 2006-05-25 07:51:28 . . . 60.2 32 

f09 52 30 2006-05-25 09:24:17 7.7±0.2 59.1 75 

-46 53 04 2006-05-25 07:51:28 . . . 60.2 32 

-36 45 42 2006-05-25 03:25:52 11.9±2.1 56.2 34 

-46 53 04 2006-08-08 08:13:56 ... 57.2 93 

-21 45 28 2006-08-08 10:09:22 10.6±1 62.9 66 



calculated from each of the 3 calibrators (see Table A.l). 
The calibrated visibilities were derived with the MIA and 
EWS packages and exhibit very good agreement between 
the data sets. The error bars on the calibrated visibilities 
are typically between 10 and 15%. The mean of the equiv- 
alent diameters of the 3 fits is 28.8 and the corresponding 
error is 0.6 mas. 




7) Syr 









I" 

a 








> 




9 10 11 1? 
Wavelength (^m) 



9 10 11 12 
Wavclcnq:h (/im) 



Fig. A. 2. Fits of equiva- 
lent uniform disks (solid 
lines) to the calibrated 
spectrally-dispersed visibil- 
ities of /3 Gruis (error bars) 
calculated from each of the 
3 calibrators : ePeg, r/Sgr, 
r 4 Eri. 



Table A. 2 lists the angular diameters of f3 Gruis found 
in the literature and those derived in this paper using the 
spectrometric and interferometric methods. 



Table A. 2. Determination of the angular diameter (in mas) of 
(5 Gruis with different methods. 



Reference 


Value (mas) 


Hcras ct al. (2002) 


28.11±2.82 


Bedding ct al. (1994) 


27±3 


Spectrometric method (this work) 


26.8±1.3 


Interferometric method (this work) 


28.8±0.6 



Diameters determined from these last two methods will 
be taken into account to reduce calibrated visibilities of 
7T 1 Gruis. The best set of calibrated visibility data will then 
be the one giving the most reliable parameter values, for 



which we will be able to simultaneously fit the source spec- 
trum and visibilities. 
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